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Abstract. The aim of this work is to prove by a perturbation method the ex- 
istence of solutions of the coupled Einstein-Dirac equations for a static, spheri- 
cally symmetric system of two fermions in a singlet spinor state. We relate the 
solutions of our equations to those of the nonlinear Choquard equation and 
we show that the nondegenerate solution of Choquard's equation generates 
solutions for Einstein-Dirac equations. 



1. Introduction 

In this paper, we study the coupled Einstein-Dirac equations for a static, spher- 
icahy symmetric system of two fermions in a singfet spinor state. Using numerical 
methods, F. Finster, J. SmoUer and ST. Yau found, in IJ, particlelike solutions; 
our goal is to give a rigorous proof of their existence by a perturbation method 0. 

The Einstein-Dirac equations take the form 



(1.2) R) ~ -R6] = -StiT} 



(1.1) {D-m)ip = Q 

1 

2' 

where D denotes the Dirac operator, tp is the wave function of a fermion of mass 
TO, i?* is the Ricci curvature tensor, R indicates the scalar curvature and, finally, 

is the energy-momentum tensor of the Dirac particle. 
In [T], Finster, SmoUer and Yau work with the Dirac operator into a static, spher- 
ically symmetric space-time where the metric, in polar coordinates (t,r,i!),(p), is 
given by 

11 



(1.3) g,j = diag 



(1.4) g'^ ^diag{T^-A,-\,- ^ 



1^ 

with A = A(r), T = T{r) positive functions; so, the Dirac operator can be written 



as 



(1.5) D^i^'dt + Y [idr + -^[i-A-^'^]-~]+i^'>d^ + i^'Pd^ 



Date: October 2, 2009. 

^After completing this work, we learned from professor Joel SmoUer that Erik J. Bird had 
proved the existence of small solutions of the Einstein-Dirac equations in his doctoral thesis in 
2005 [2]. His method is quite different from ours: he uses Schauder's fixed point theorem. 
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with 




(1.6) 




(1.7) 


Y 


(1.8) 




(1.9) 


7^ 



\/A {j^ cos I? + 7^ sin d cos + ^ sin sin (/?) 
- (—7^ sin i9 + 7^ cos -d cos + 7'^ cos -d sin (p) 



-7^ sin + 7^ cos 93) 



r sin d 

where 7* are the Dirac matrices in Minkowski space (see [T]). 

Moreover, Finster, SmoUer and Yau are looking for solutions taking the form 



(1.10) 



$1 



\ 



J 



where cr'" — {a^ cos "d + d'^ sin d cos ip + sin sin 1^9) is a linear combination of the 
Pauli matrices ct* and <i>i(r), $2(?') are radial real functions. 

We remind also that the energy-momentum tensor is obtained as the variation 
of the classical Dirac action 



S ■ 



ip{D — m)%l)\/\g\ d x 



and takes the form 
1 



diag 2wr2|4>p, ~2ujT'^\<^\'^ + 4T-$i$2 + SmT ($? - $i) 



1 1 

-2r-$i$2,-2T-$i$2 
r r 



(see [T] for more details). 

In this case, the coupled Einstein-Dirac equations can be written as 



(1.11) 


^/A<^>[ 


= i$i - (tjT + to)$2 
r 




(1.12) 


\/l$'2 


(uT - m)$i - i<l>2 
r 




(1.13) 


rA' 


= \- A~ leiTLuT'^ ($2 - 






T' 

2rA— 
T 


= A-1- IGttluT^ ($2 - 


^$2) + 327riT-I>i<I>2 
r 


(1.14) 




+167rmr(<I>^ - <^l) 





with the normalization condition 
(1.15) 



T 1 
|$|2J=dr = -i. 
^/A 47r 



In order that the metric be asymptotically Minkowskian, Finster, SmoUcr and Yau 
assume that 

lim T{r) = 1. 

r — >QO 

Finally, they also require that the solutions have finite (ADM) mass; namely 

r 

lim -(1 - A{r)) < 00. 

r— >oo 2 



PERTURBATION METHOD FOR EINSTEIN-DIRAC EQUATIONS 



3 



In this paper, we will prove the existence of solutions of (|l.im.2p in the form 
()1.10|) by a perturbation method. 

In particular, we follow the idea described by Ounaies in [5] (see also ^ for a 
rigorous existence proof of nonlinear Dirac solitons based on Ounaies' approach). 
Ounaies, by a perturbation parameter, relates the solutions of a nonlinear Dirac 
equation to those of nonlinear Schrodinger equation. Imitating the idea of Ounaies, 
we relate the solutions of ours equations to those of nonlinear Choquard's equation 
(see [3], [5] for more details on Choquard's equation) and we obtain the following 
result. 

Theorem 1.1. Given < tu < m such that m — uj is sujjiciently small, there exists 
a non trivial solution of II 1 . ll'pTT^ . 

In Section [2l we solve the Einstein-Dirac equations by means of the perturbation 
method suggested by Ounaies; in particular in the first subsection we describe 
a useful rescaling and some properties of the operators involved, whereas in the 
second subsection we prove the existence of solutions generated by the solution of 
the Choquard equation. 



2. Perturbation method for Einstein-Dirac equations 

First of all, we observe that writing T(r) = 1 + t{r) and using equation (|1.13p . 
the coupled Einstein-Dirac equations become 

(2.1) VA^\ = i$i - (cj + m)$2 - ^i$2 

r 

(2.2) \/]4$2 = (w - m)$i + - 

r 

2rAt' = (A- 1)(1 + <) - 167rw(l + i)^ ($? + $2) + 

(2.3) +327ri(l + t)2$i$2 + 167rm(l + tf {<^>l - ^j) 

where 

(2.4) A{r) ^l-^J\l + tis)r ($i(5)2 + a>2(.)2) ds:^l- ^Q{r). 

Furthermore, because we want A{r) > 0, we have that the following condition must 
be satisfied 

(«) — <Ti^ 

r lOTTOJ 

for all r G (0, 00). 

Now, to find a solution of the equations (|2.HI2.3p . we exploit the idea used by 
Ounaies in In particular, we proceed as follow: in a first step we use a rescaling 
argument to transform (|2.m2.3p in a perturbed system of the form 



(2-6) \ ^ A (e, ^, X, t) + + <P - "i^t + K2 (e, X, t) = 
A (e, ^, X, r) f r + «^ /; ^2 d5 + (e, x, r) = 
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where (p, x, r : (0, oo) M. 

Second, we relate the sohitions of (j2.6p to those of the nonUnear system 

(2.7) <! X(r)-5k(^^-f^) 

r(r) = 87rm r„°° — — r ds. 

^ / JO max(r,s) 



ds ) iy9 = 



We remark that is a solution of ()2.7p if and only if u{x) 
nonlinear Choquard equation 

(/«3 ^ dy 



solves the 



(2.8) 



-Am + 2mu - 



in 



\x-v\ 

To prove this fact it's enough to remind that for a radial function p, 

^2 ^j. 



Ap 



2 

r —p 

ar 



and 



(I-Kp) (a;) ===4^ 



max(r, s) 



ds 



with r = |a;|. 

We observe also that if we write 



u{x) 
v{x) 



( 



\ 



with r : 
(2.9) 



V 



(p(r) 
zx(r)fT''' 

\x\, (iy9,x) is a solution of (|2.7p if and only if (m(x), i;(a;)) solve 
-Au + 2mu - 4m3 /jj3 Hlil^ dy) w = w = 



in M3 where ctV = Ya=i ^'di. 

It's well known that Choquard's equation (|2.8p has a unique radial, positive solution 
ug with J |moP = N for some iV > given. Furthermore, Uq is infinitely differen- 
tiable and goes to zero at infinity; more precisely there exist some positive constants 
C5,r, such that \D'^ (uo)| < Q.r, exp(-(5|a;|) for x e R? . At last, mq G H'^{R?) is a ra- 
dial nondegeneratc solution; by this we mean that the linearization of (|2.8p around 
uo has a trivial nuUspace in L^(R'^). In particular, the linear operator C given by 

|2 \ 



£^ = -AC + 2m^-4m^ 




uo{y)Y 
x-y\ 



i{y)uo{y) 



y\ 



dy uo 



see [4], [5], [3] for more details), 
are the zeros of a operator D : 



satisfies ker C — {0} when C is restricted to L^( 
The main idea is that the solutions of (|2.6p 
M X X X Xr ^ X X Yt- If we denote by D^^^,r{£, 'P, Xi ''') the derivative 
of D{e, ■,-,■), by {(poiXo^To) the ground state solution of (|2.7p and we observe 
that D^p^^^ri^i 'fo, Xo, ''0) is an isomorphism, the application of the implicit function 
theorem (see [7]) yields the following result, which is equivalent to theorem ll.il 



Theorem 2.1. Let {(Pq,xo,tq) he the ground state solution of Ii2.7\ ), then there 
exists S > and a function rj G C((0, S),Xtp x X^ x Xt) such that 77(0) = {(po, xo, tq) 
and (e, 77(e)) is a solution of 112. 6\) . for < e < S. 
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2.1. Rescaling. In this subsection we are going to introduce the new variable 
X,t) such that <i>i(r) — a(^(Ar), $2(?') — Pxi^f) a-nd t{r) = 7r(Ar), where 

$1, $2, i satisfy (|2.1j|2.3p and a, /?, 7, A > are constants to be chosen later. 
Using the explicit expressions of A, given in (|2.4p . we have 



(2.10) := A^^p,^{ip,x,T) 



It's now clear that if $1, $2, t satisfy ()2.m2.3|) . then (p, x, t satisfy the system 



r ./A^^fj-^p-flp+{m + uj)x + u:^Tx^Q 
(2.11) I ,/A:j^^f^x+lx+fxi^~^W-jk^rp = ^ 



with 



i^a,/3,7(^,X,T) = ^^ n (l + 7r)2 L^^ (^^x) j ds\ (l+7r) + 



+167rc.^(l + 7T)3 _32^Ml(i+^^)2^^ + 

^2 / 2 ' 

— loTrm — (1 + 7r) — —X 

7 \ 



By adding the conditions -^{m — oj) = 1, ^ = 1, ^ = 1, ^ = 1 and m — uj > 0, 

we obtain a — (m — cjY^^ , A = (m — w)"'^/^, (3 — m — uj and ^ = m ^ uj. 
Denoting e = m — tj, (|2.1ip is equivalent to 



\/A {e,ip,x,T)£ip - lip + 2mx + Ki (e,(^,x,T) = 
(2.12) <( (e. X,t)£x + Ix + 'P- m^r + {e, (p, x, r) = 

A{e,ip,x,T)£T + ^j;;^^ds + Ks{e,ip,x,r)^0 



where A (e, x, r) , if 1 (e, (/3, x, t) , if2 (e, V, X> t) and K3 (e, x, t) are defined by 

(2.13) A {e, p,x,r)^l- ^^"^"^ ~ (1 + er)^ (^^ + ex') ds; 

^ Jo 

(2.14) Ki (e, X, r) = -ex + £("i - £)tx; 

(2.15) if2 (e,'^,X,T) = ETV'; 
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and 

K3{e,'f,x,r) = X^ds+—^ t {^^^ + ex^) ds + 

+— ^ / r^^^+ex^) ds + 
^ Jo 

-'-^{[i^ + ^rfi^^+ex') d.)r + 

+167rme— + 87rme(3T + Ser^ + eV^)^^ ±-L + 

r r 

_8;re(l + err^^^±f^ - I6ne{l + er)'^ + 

(2.16) -87rme(2T + er2)-!^-— 

r 

For £ = 0, p.l2p becomes 

|.<^ _ + 2mx - 




dr ' 



4zT + ^ /; ds = 



+ 2m^ - 16 W ^^^^^^ ds) ^ = 

Then, we denote by (930, XO: ''0) a solution of (|2.18p : in particular 

^^^^sj ^ r_d_ m\ 

2rn dr \ r J 

f°° m"^ 
To(r') = Svrm / -ds. 



Iq max(r, s) 

Now, to obtain a solution of (I2.12p from {ipQ, xq, tq), we define the operators Li : 
RxX^xX^xXr Y^, L2 : RxX^xX^xX^ Y^, L3 : RxX^xX^xX^ Y^ 
and D : R X X^ X X^ X Xr ^ Y^ X Y^ X Yr by 

-Id X 1 

)- — <y9 + 2m h 

r dr r-^ r r 

-1 . X , f f . ^ 



Li{e,(p,x,T) = \/A {e,(p,x,T)-- — (p + 2m-- + -Ki {s,(p,x,T) 



L2{£,V-,X,t) = ^/A{e^^p~xy)--rX + m-r + -K2{e,tp,x,T) 

r dr j.^ J- J. J. 

T I \ A I \ d Bttto r 2 , , \ 

L3(£,'^,X,T) = A{e,(p,x,T) —T + J if ds + K3[e,ip,x,T) 

and 

Die,(p,x,T) = {Li{e,(p,x,T),L2{e,(p,x,T),L3{e,ip,x,T)) 
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where 



= |(^:(0,oo)-. 

= |x:(0,oo)^ 

X, = |t:(0,oo)^ 

K = 



y(l-^i) ( 1 

\x\ \ 
X(l^l) 



lim r(r) 0, -^r € ^^((0, oo), dr) i 
r^oo dr J 



L1((0,oo), (ir). 
Furthermore we define the following norms: 

<pi\x\) 



llxllx, 

\\r\\x^ 



X{\x\) ^r ( 1 



^1(113) 



Li((0,oo),dr) 



It's well known that 



L°°((0,oo),rfr). 



Moreover, using Hardy's inequality 

|2 



Jw' \x\ Jr3 



we get the following properties: 

pG ((0, oo), dr) ^ ((0, oo), dr) 
^eL2((0,oo),dr). 



(2.19) 



ype x^yp£ x^. 

Since the operator ^(e, (/j,x, t) must be strictly positive, we consider B^, B-^, 
Bt, defined as the balls of the spaces X^,Xy^,Xr, and ei,e2, depending on m and 
on the radius of B^, B^,Br, such that 



167r(m-e)£ /"",_, / 2 2\ 

1 ^ ^ / 1 + erf ^2 ^ 



> (5 > 



for all (e, (p, x, t) € (—£1, £2) x x x i?^. The existence of £1, £2 is assured by 
the fact that tp, x, t are bounded; in particular, if £ > 0, 

167r(m — £)£ 



1 - 
> 1 



■/ {l + eTr{^^+ex')ds> 
Jo 



+ 



20me - 8me' (5 ||r||^^ \M%^ + ||x||^J 

4ms' \\t\\^^ (5 llrll;,^ + 4 \\x\\l^) - 8me' ||r||^^ ||x||^^ 
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then there exists £2 > such that A{e, ip, Xj ''') > for all e € [0, £2). In the same 
way, if £ < 0, 



1 - 



167r(TO — £)£ 



(1 + £t)2 (<^2 g^2^ > 



> 1 - 8me' WxWl^ - 8m|£|3 (l + 2 ||t||^J + 

—8me^ llxl'^ 



(2 + l||r||^J 



8Me\'\\r\\lj\xfx^, 

then there exists £1 > such that A{e, (p, Xi ''') > for all e E (—£1, 0). 

Lemma 2.2. T/ie operators Li,L2 E ((— £i,£2) x x x Bt-,Y^p) and L3 E 
Ci((-ei,£2) X B^x B^x Br,Yr). 

Before starting the proof of the lemma we observe that for a radial function p 
such that ^ E Hi (W^) we have 



(2.20) \p{r)\ < 

We remind that (M.^) = {u E 



d f pir) 



dr \ r 



u IS r, 



adial } . 



Proof. We begin with L3; first, we have to prove that it is well defined in Yr 
L^((0, 00), dr). Wc remark that 



\L3{£,'P,X,t)\ < Ci 



dr 



% fk^+ex^l ds + ^W^+ex^\ 



C, 



where Ci, C2, C3, C4, C5 are positive constants and, by definition, we have that 

freLi((0,(X)),dr). 

Next, we have 



— / \'p^ + ex^ I dsdr — 



°° Iv^'+ex'! 



ds < +00, 



/o 'Jo Jo 
using Holder's inequality, then |(^^ + £x^| ds E Y^. In the same way, we can 

conclude that i {ip^ + £x^) , H'^^ " s^^) ^ 
Finally, 







X 




r 


L2((0,oo)) 


r 


L2((o,oo)) 



< +00 



thanks to fU^ . then ^ipx eY^-. 

Now, we have to prove that (£, 93, x, t) is C^; by classical arguments, it's enough 
to show that for {hi, /i2, h^) E B^p x B^ x Br 

d 

■q^ {L3 {e,ip,x,'r))hi E Yr, 
d 

— {L3 {e,ip,x,'r))h2 E Yr, 
ox 



^ (-^3 (e,¥',X,T)) E Yr. 
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We begin with ^ (L3 (e, ip, X, r)), 
d 



dip 



d 



— {L3 (e, <p, X, r)) hi= (—{A (e, <p, x, r)) hi ) — r + 



dtp 



167r(m — e) 



Jo 



(phi ds] {1 + £t) + 



+167r(m - + er)^^ - 167r£(l + srf^ + 



-167rm(l + st) 



2</?/ll 



for S^{Lz {e,ip,x,T)), 



d 



{L3 {£,<p,X,r))h2 



dx 

167r(m — e)e 



d 



{A{s,ip,x,T))h2 ) — T + 



dx 

(1 + eT)^xh2 ds]{l+ et) 



+ 



+167r(m - e)e{l + srf^ - 16ws{l + srf^ + 

-167rme(l + eT)2^ 
r 



and, finally, 



d 



d 



— (L3 (e, X, t)) /is = ^ (e, X, r)) h3]—T + 



dr 



d 



dr 



d 



+ 



+ 



16w{m — e)e 



87r(m — e)e 



^ (1 + £r)/i3 {'p'' + £X') dsj (1 + £t) + 
{1 + eT)^p^ + ex^) ds] /is + 



+247r(m - £)e(l + er)^^'^^ ^ ^^^^ ^ 



/i3 - 327re^(l + £T)^/i3 + 

r 



-167rme(l + er)^ ^^/i3- 

r 

First of all, we remark that if (p,hi € B^, Xih^ G -B^ and t, ^,3 € Bt, then 
^ (^(£,<P,X,t))/ii, (A(£,(^,x,t))/i2 and ^ (A(e, v?,x,t)) /is are bounded. 
So, we have that 



dp 



dx 
dr 



< Ci 

< C5 

< Cg 



dr 
d_^ 
dr 
d ^ 
dr 



+ 



.21 ds]+Cr^-^ + Cs 



\<ph2\ 



13- 



\v>x\ 



+ c 



-/is 

dr 



12 21 



r 

„2 I „,,2| 



+ 



14- 
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with d positive constants. With exactly the same arguments used above, we 
conclude that 



d_ 

d_ 

dx 
d_ 



{L3 {£,V,X,r))hi 
{L3 {e,'P,X,T))h2 
{L3 {e,'P,X,T))h3 



dr < +00 
dr < +00 
dr < +00 



if {e,(p,x,T) G (-61,82) X X X Br, and (/ii,/i2,ft3) G B^ x B^ x Br 



Furthermore 



dLs dL. 



, and are continuous; thus the proof for L3. 



We consider now Li, first, we have to prove that it is well defined in Y^. We 
observe that 



\Li (e,</?,x,r)| < Ci 



1 d 
-—if 
r dr 



C2 



with Ci,C2 positive constants then, Li {e,ip,x,T) £ (R^), thanks to conditions 

Now, we have to prove that Li (e, ip, x, t) is C^; by classical arguments, it's enough 
to show that for {hi, /12, h^) £ B^p x By, x Br 



d_ 

dip 
d_ 

dx 
d_ 
dr 



{Li {e,ip,x,T)) hi e Y^, 
[Li {e,ip,x,T))h2 G Y^, 
{Li {e,(p,x,T))h3 G Y^. 



By a straightforward computation, we find out 

dip ^ 2 \dp ^ ) r dr^ r dr ^ 



h2 

r 



dLi , 1,1/9/ dA , \ 1 d , 

1 ,_i/2 fdA^ \\d , , V 



h2 

r 



and, using the positivity of A, 



dLi^ 
dip 


< 




1 d 
-—tp 
r dr 


+ C2 


dLi 
ox 


< 


C3 


1 d 
-—tp 
r dr 


+ C4 


dLs, 


< 


C5 


I d 
-—p 

r dr 


+ Ce 



1 d 

r dr 

r 

X 
r 



hi 
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with Ci positive constants. Then, we can conclude that 

2 

dx < +00 





d 


/ 






d 


/ 


dx 




d 


/ 


97 



(^3 {£.^,X,T))hi 
(Li (£,V,X,T))h2 



dx < +00 



dx < 



if (e, (/?, x,t) e (-61,62) X X X Br, and (fti, /12, G B^ x B^ x Br. 
Furthermore and are continuous; thus the proof for Li and with the 

same arguments for L2- 

□ 

2.2. Branches generated by solutions of Choquard equation. In this sub- 
section, we show that a solution (/)o = ((ySo, Xo, 'Ha) of (|2.7p can generate a local 
branch of solutions of (|2.6p . 

First, we linearize the operator D on r) around (0,(/)o) 

/ i:^/i-^ + 2m^ \ 



1 _d_ 

1 ''-fc + 4 



V 



Lemma 2.3. We define the operator V : x x Y^, by 



r ar < < j. \j 

^l + ^Io^ohds 

Now, if we prove that D^^^^r{0, (po) is an isomorphism, the implicit function theorem 
can be applied and we can find solutions of (|2.6p near the ground state ^o- 



r dr^ 



X 



then V is an isomorphism of X^ x X^ onto Y^ x Y^. 

This lemma is obvious if we remind that L^(M'', C*) can be written as the direct 
sum of partial wave subspaces and that the Dirac operator leaves invariant all these 
subspaces (see [8J). So, thanks to lemma 2.3 of [B], we know that V : 
(R3^C2) ^ L2 (R3^c2^ X (-r3^c2) defined by 



V{u,v) 



iaVu + 2mv 
—ia'Vv + u 



is an isomorphism of (R^, C^) x (R^ C^) onto (M^ C^) x (M?, C^) and 
then V is an isomorphism of each partial wave subspace. In particular, V coincide 
with V on the partial wave subspace X^p x X^ . 

Lemma 2.4. We define the operator W : X^ x X^ x Xr Y^ x Y^ x Yr, by 

T dr T 



r dr 



k 



\ 



-2-/ 

dr 



W{h,k,l) 

then W is an isomorphism of X^ x X^ x X^ onto Y^p x Y^ x Yr- 



J 
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Proof. First we prove that W is one to one. We observe that W{h, k,l) = if and 
only if {h, k, I) satisfies 

r dr r 

-T-fc + -^ + --m^Z = 

r dr > > J- J- 

dr 

in X X Yr- In particular, we must have I = and (/i, k) solution of 

--fh- \ + 2m^ = 

r dr r 

-4-k+\ + - = 

r dr r^ r 

that is equivalent to V{h, k) — 0. So, thank to lemma [^751 h = k = and ly is one 
to one in Y^ x Y^ xY^. 

Secondly, we have to prove that for / = (/i, /2, /s) &Y^xY^xYr, there exists 
{h, k, I) G X X Xr such that W{h, k, I) = f. This means that the system 

-T-^ - 4 + 2m^ = A 

r dr r -'- 

. dr^ ~ 

has a solution in x X^ x Xr for all (/i,/2,/3) G x x IV- We observe 
that V/a e i^((0,oo),dr) there exist l*{r) = - f^ds such that = /s; 

furthermore /* G X,-. So, we have to show that 

(2.21) <^ "'''^ 

has a solution in X^ x X^ for all (/i, A) G ^x- 

Now, we remark that ^l* G (M^) and then, thanks to lemma [231 has a 

solution in X^ x X^ for all (/i, A) G Y"^ x Y^. 

In conclusion W is an isomorphism of X^ x X^ x Xr onto Y"^ x K,, x F,-. 

□ 



Finally, we observe that _D^^^.t-(0, 4'o){h, fc, I) can be written as 

(2.22) A, 
with 

/ \ 

(2.23) S(h) = 



'^,x,r(0, M{h, fc, /) = M/(/i, fc, Z) + 5(/i) 



Theorem 2.5. Let <j)Q he the ground state solution of \2. 7^ , then there exists 5 > 
and a function rj G C((0, S), X^ x X^ x X^) such that 77(0) = 0o ^md D{e, rj{e)) — 
forO<s<S. 

Proof. Since D{0, (po) = and D is continuously differentiable in a neighborhood of 
{0,(j)o), to apply the implicit function theorem we have to prove that D^^^^r{0,4'o) 
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is an isomorphism of x x Xr onto x Y,. x y,-. 

We observe that D^^^^r{Q^ Mih, k,l) — if and only if {h, k, I) satisfies 

(2.24) < £k+^ + h-mhTo-mipQl^O 

that means 
(2.25) 

_ 



h + 2mh - Wnm^ ( r„°° , ds) h - 32Trm^ ( '^f , ds) ipo = 

\ JO max(r,s) J \ JO max(r,s) y ^ u 

^h- ^h + 2mk^0 

dr r 

I = 167rm '^f , ds 

Jo inax(r,s) 

Now, if we write ^{x) — -^^^jjp and we remind that (/3o(|a;|) = |x| 1*0(2;) with uq 
solution of (|2.8p . we have that {h, k, I) is a solution of (|2.25p if 



Hr) ^ 
y jA:(r)cr'' q y/ / 



satisfies 

(2.26) < 

I /• - -'g^? 

2m 

and 

(2.27) Z(.) = 4m/ dy. 



-AC + 2K - 4^3 (/^3 dj;) e _ (/«3 dy) u, = 

C 



It's well known that the unique solution of the first equation of (|2.26p in H^{M.^) is 
^ = (see [3] for more details) and that implies (^ = 1 = 0. So the unique solution 
of (|2.24p ish = k = l = and D^^^^t{0, (j>o) is one to one in X^p x X^ x X^-- 

Next, if we show that S{h) is a compact operator, we have that D^^^.t-(0, 0o) 
is a one to one operator that can be written as a sum of an isomorphism and a 
compact operator and then it's an isomorphism. 

First, we can easily see that T{h) = ^ [J^i^^hds^ is a compact operator from 
X^f, on Yt-\ in particular, we use the fact that {M?^ is compactly embedded in 
L'^ (M"^), for 2 < (7 < 6, to prove that for any bounded sequence C X^, the 

sequence {r(/i„)} C contains a Cauchy subsequence. 

Second, we have to show that the operator -^ro from X^^ to L? (M?^ is compact. If 
{^} is a bounded sequence in iJ^(M3) then {^tq} is precompact on Lf^^{^^), 
thanks to compact Sobolev embedding and, since ro(r) — > when r — > +cxd, we can 
conclude that {^tq} is precompact on L^(R3). 

So S{h) is a compact operator from X^ on x x and D^^^^ri^^M is an 
isomorphism of X^ x X^ x Xt onto Y^p x Y^ x Y.,- . 

In conclusion, we can apply the implicit function theorem to find that there 
exists 5 > and a function 77 e C((0, (5),X(^ x X-^ x Xr) such that 77(0) = (/>o and 
D{s,r){s)) = for < e < (5. 

□ 
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